Long-range-correlated disorder in graphene 
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We study transport of two-dimensional quasi-relativistic electronic excitations in graphene in the 
presence of static long-range-correlated random scalar and vector potentials. Using a combination of 
perturbation theory and path-integral techniques, we estimate scattering rates which control Drude 
conductivity, magneto-transport, and Friedel oscillations in the ballistic regime of large quasiparticle 
energies. We also discuss properties of zero-energy states and pertinent localization scenarios. 



In recent years, the theory of electron localization in 
two dimensions (2D) has been extended to the situations 
where, instead of (or in addition to) a random poten- 
tial (RP), there exists a random magnetic field (RMF). 
This type of problems emerge in the context of compress- 
ible Quantum Hall states 1 , finite-temperature dynamics 
of spin liquid states in Mott insulators 2 -, and vortex line 
liquid phase in high-T c cuprates 3 -, to name a few. 

For example, in the latter case fermionic excitations 
possess quasi-relativistic (albeit not necessarily rotation- 
ally invariant) kinematics, and the roles of RP and RMF 
are played by the Doppler shift due to a circulating su- 
percurrent and concomitant Berry phase, respectively 3 . 

Another relevant application of this disorder model can 
be found in graphene which has received a lot of attention 
lately. Although a number of different (anti)localization- 
related phenomena have already been discussed^, these 
studies were largely limited to the case of short-ranged 
(albeit, possibly, arbitrarily strong) disorder. 

Recently, the effect of long-range Coulomb impurities 
residing in the Si02 substrate was invoked^ to explain 
the experimentally observed linear dependence of the 
conductivity on electron density^. 

It has also been pointed out that in graphene, be- 
sides the Coulomb ("scalar") RP, there exists a "vec- 
tor" (RMF) disorder representing the effect of disclina- 
tions (isolated pentagon- and heptagon-rings), disloca- 
tions (pairs of adjacent pentagons and heptagons), and 
Stone- Wales defects (double pairs)£. The presence of 
such structural defects in free-standing graphene is in- 
evitable due to the intrinsic thermodynamic instability 
of 2D crystals. 

Nevertheless, in most of the previous work this RMF 
was repeatedly treated as being short-range in terms of 
the vector potential (rather than the magnetic field it- 
self), which simplification facilitates the use of a powerful 
machinery of renormalization group 8 and 2D conformal 
field theory^. 

In the present Letter, we compare the effects of long- 
range-correlated RP and RMF on electron transport in 
graphene and contrast the results with those pertaining 
to the conventional 2DEG with parabolic electron dis- 
persion. 

The Dirac-like quasiparticles near two conical points 
(K and K') in the hexagonal Brillouin zone of graphene 



cn be described in terms of the (retarded) Green function 



(e + S fl )7 + tryp 
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where v is the Fermi velocity and the 4x4 7-matrices 
7 M = [il (T 3) 1 (g) a 2 , — (T3 <g> <7i) act in the space of the 
Dirac bi-spinors ip = (ip K (A),tpK(B),ipK'{A),ipK'{B)) 
composed of the values of the electron wave function on 
the A and B sublattices of the bipartite hexagonal lattice 
of graphene. 

The Dirac fermions are subject to random scalar (s) 
and vector (v) fields whose spatial correlations are con- 
trolled by the Gaussian averages 

(a fl (q)a v (-q)) = S f _ l0 S l/0 w s (q)+S f _ ll S l , j w v (q) fs^ - 

q (2) 

The variances w s , v (q) = 47r 2 r S!t ,/q 2 are proportional to 
the arcal densities of the Coulomb impurities (T s = g 2 rii 
where g — e 2 / £qv is the Coulomb interaction parameter) 
and topological defects (T v ~ n<j), respectively. 

Strictly speaking, in the case of graphene the vector 
field a becomes a 4 x 4-matrix proportional to (1 <g> 03)^, 
thus endowing Eq.(2) with the structure of a 16 x 16- 
matrix. It appears, however, that this technical compli- 
cation does not affect any of the results presented below 
and, therefore, we will use Eq.(2) as is. 

A preliminary insight into the problem can be gained 
by attempting to compute a fermion self-energy in the 
framework of the customary self-consistent Born approx- 
imation (SCBA) 
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In the case of RP, the singular behavior of iv s (q) at 
q — > gets replaced by w s (q) = T s /(q + k) 2 due to 
the Debye screening. The corresponding Debye momen- 
tum k = 4:gpF proportional to the Fermi surface radius 
p F = (true) 1 / 2 , where n e is the density of excess elec- 
trons with respect to half-filling, renders the correspond- 
ing quasiparticle width finite 



7 S = ImTrj £K(e,e/v) 
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By contrast, a calculation of the self-energy associated 
with RMF is impeded by the lack of screening for the 
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random vector potential which results in an infrared di- 
vergence 



lv = ImTr%± R (e, e/v) ~ vT^^ 



(5) 



where L is the size of the system. 

This intrinsic divergence can not be avoided even if 
one proceeds beyond the SCBA, for it stems from the 
non-gauge invariant nature of the fermion Green function 
G R (t, r) for r 7^ 0. Taken at its face value, the divergent 
self-energy (5) is indicative of a strongly non-Lorentzian 
form of the RMF-averaged Green functionli. 

It is worth mentioning that, unlike in the previ- 
ously discussed examples of the RMF problem involv- 
ing some auxiliary fermions which are different from 
physical electron*!*^, in the present case there is no 
physical ground for replacing the original Green func- 
tion (1) with a gauge-invariant amplitude such as, e.g., 
G R (t,r) =< ^(i,r)exp(-i J c a(r')dr')^(0,0) >. 

Nonetheless, Eq.(3) can be used to evaluate the Drudc 
transport rates in both, RP and RMF, cases. Inserting 
the factor 1 — cos 9 related to the transferred momentum 
q = 2p sin 9/2 into the integrand in Eq.(3), and putting 
e = ep = vpp, one obtains the first-order Born estimates 



S(e F - v\p + q\)w a {q) sin 2 9 (6) 



ponent. In contrast to various approximate (e.g., Bloch- 
Nordsieck) representations, the integration over the mo- 
mentum p(r) conjugate to the spatial coordinate r(r) al- 
lows one to account for the spinor structure of the fermion 
propagator (1) exactly. 

Averaging over the disorder variables introduces a 
product of the " Debye- Waller" attenuation factors into 
the integrand in Eq.(8) (here = (l,dr/dr)) 

u„(T)u„(T'){o„(q)a„(-q)) e "><'<'>-'l''»] (10) 

Despite its generally gauge-non-invariant nature, the am- 
plitude (8) does appear to be gauge-invariant for closed 
trajectories (r = 0), thus allowing one to analyze such 
magneto-transport effects as dHvA (SdH) oscillations of 
the density of states (hence, magnetization, etc.) and 
conductivity in a weak uniform external magnetic field 

£<r S! „. 

Evaluating Eq.(8) on the semiclassical trajectories 
which dominate the path integral and correspond to mul- 
tiple repetitions of the Larmor orbit of radius R c = e/vB, 
one obtains a field-dependent density of states 



u(e\B)=u(e\0) 



a 2KinA(e)-n 2 8S 1 (e) 



(11) 



which both turn out to be finite 

Is min[l, —t\, 

e F g l 
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(7) 



1/2 

and inversely proportional to n e , thereby giving rise 
to the Drude conductivity a oc n e , in agreement with 
experiment^ (note that the experimentally relevant value 
of the Coulomb coupling is g ~ 1) and in contrast to the 
situation in the conventional 2DEG where ^ ~ T a /m 

1/2 

contain the band mass m, thus yielding a oc n e ■ 

To make a further progress, we employ a path-integral 
representation of the Dirac fermion Green function which 
was devised in Ref and has been previously applied to 
the problem of quasiparticle transport in the vortex line 
liquid phase of the cupratesi^ 



/>oo pr(t)—r 

G fl (e,rMr)) = dt Dr Dp 

JO Jr(0)=0 

/■* dr 
exp[iS (t) + i / dT(a (r) - — (A + a(r))] (8) 

Jo dT 



where represents external field (if any) and the free 
Dirac action reads 



S (t) = rfr[e 7o + p(— -7)] 



(9) 



In Eq.(8), the usual ordering of 7-matrices with respect to 
the proper time r must be performed according to the or- 
der of their appearance in the series expansion of the ex- 



where A(e) = ire 2 /B is the area of the Larmor orbit and 
the attenuation factor contributes as 

^i(e) = - £ lnW a =n[^ + ^} (12) 



which allows one to identify the proper cyclotron rates 
7^~< /2 , T^-feVr.) 1 / 4 (13) 

associated with the linear and quadratic (in powers of 
1/B) Dingle plots, respectively. 

These results should be contrasted with their "non- 
relativistic" counterparts (7^ yc/ ~ T s /m, j^ ycl ~ 
{(Ty/m) 1 / 2 ). In combination with the transport rates 
(7), distinct energy (hence, electron density) dependences 
of the rates (13) make it possible not only to distinguish 
between the conventional 2DEG and graphene, but also 
to discriminate between the RP and RMF mechanisms 
of elastic scattering (experimentally measured rates (13) 
would be naturally estimated for e = ep). 

With the use of the Poisson formula, Eq.(ll) can be 
cast in the form 



v{e\B) oc ^]exp[ 
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v 2 B{^ vcl ) 2 + (lv yc Y 



(14) 



showing an inhomogeneous broadening of the relativistic 
Landau levels positioned at ui n = ±v(2Bn) 1 ^ 2 . 

The path integral technique is also well suited for an- 
alyzing various gauge-invariant two-particle amplitudes. 
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For one, the ballistic Drude conductivity manifesting the 
transport rates (7) can be found by computing the av- 
erage < G R (e, r)G (e, — r) > in the semiclassical ap- 
proximation. Expanding the corresponding path inte- 
gral about a proper semiclassical trajectory (see below), 
one can systematically study ballistic corrections to the 
Drude conductivity 14 . Due to the long-range nature of 
the RMF disorder, these corrections do not appear to be 
logarithmic and, therefore, are not readily amenable to a 
simple resummation via renormalization group (cf. with 
the case of short-range disorder—) . 

Another example of a relevant two-particle amplitude 
is given by the correlation function of the electron wave- 
functions' amplitudes related to the average^ 

/•oc />i\(ti)=±r 

< G R (e,r)G R (e,-r) >= TT / dU Dr z D Pl 

8=1,2^° J ri (o)=o 
e iSo( tl ) e iSo(fe) "Q Y[ W a {n(n)-rjfa)) (15) 

a— s,d i.j — 1,2 

Observe that the product of the W^-factors yields an ex- 
ponent of the Amperian area of a contour formed by the 
trajectories ri(r) and — r 2 (r). 

In the ballistic regime, the integral (15) receives its 
main contribution from pairs of trajectories with single- 
valued projections onto the semiclassical straight-path 
trajectory r«(T) = rr/t. 

Proceeding by analogy with the earlier calculations in 
the case of the conventional 2DEG^^, we separate out 
the coordinate variables onto the "center of mass" and 
relative ones (r^ = ri ± r2, — pi ± P2), expand up 
to the second order in r and p , and integrate over 
all the variables except for the transverse deviation from 
the straight-path a;j(r), thus arriving at the disorder- 
induced correction to the free Dirac action 



can be presented in the form 



5s (o,oh 

g v (0,Q\v) = 
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(18) 



In the ballistic regime, the scaling functions f s ,v(z) ap- 
proach the parabolic cylinder and Airy functions, respec- 
tively. 

Computing a real-space asymptotic behavior of the ir- 
reducible part of the wavefunction amplitudes' correlator 



L 4 < |^ 2 (r)^ 2 (0)| > -1 



< ImG R (e, r)ImG R (e, -r) > 



(^) 1/2 cos(2er)e- I " ) '° (19) 



one can identify the rates which control the decay of its 
Friedel-type oscillations 



SS 2 (t) = f dr\T, 

Jo 



e(xl(r)) 2 + T v v\xl(r)\] (16) 



Comparing typical values of the free fermion action and 
the correction (16), we find that the condition So 3> SS, 
under which the path integral (15) would be dominated 
by the trajectories close to r^(r), is readily satisfied in 

1 /2 

the ballistic regime (e ^> vTd )■ 

As in the previous wor k 2 i 13 i 15 , it can be shown that the 
amplitude (15) is given by a product of the free Dirac 
propagator with ImG R (e,r) ss e(l + 7 f)Jo( er ') and a 
Fourier transfrom (over the variable e — p) of the Green 
function g(x, x'\e,p)\ x=X ' = o of the ID equation of motion 
in the direction perpendicular to the classical trajectory 

[v 2 dl+(e+ivT v \x\+ieT s x 2 ) 2 -v 2 p 2 }g(x,x'\e,p) = 5(x-x') 

The imaginary effective potential appearing in (17) re- 
strains the transverse Dirac fermion's motion, unlike a 
real potential which allows for the Klein tunneling. 
In the pure RP or RMF case, the solution of Eq.(17) 



r 2/3 
J- V 



(20) 



and contrast them against their non-relativistic counter- 
parts ( 7 f ~ 7 f and 7 ™ ~ r 2 /^ 1 / 3 /m 2 / 3 ). 

These rates (equivalently, length scales) can also be 
manifested by other types of Friedel oscillations, such as 
that of the electron density profile induced by an iso- 
lated impurity or the RKKY interaction between a pair 
of magnetic ions. Namely, at distances r ~ »/)™ the 
previously found behavior Sp(r) cx cos(2er)/r 3 — changes 
to 



Sp(r) oc (2iL_) 1 / 2 cos(2er)e- r ^ 



(21) 



1/2 

In the complementary low energy limit (e < vT a ), 
the above eikonal-type approach ceases to be applicable. 
Nonetheless, one can still gain some insight into the local- 
ization properties of the system in question by focusing 
on zero-energy states (if any). 

In the (apparently, more challenging) case of a pure 
RMF, these states can be explicitly constructed in the 
form 



e 0W 

■0±(r) (X (1± 70 ) ( o _0( r ) 



(22) 



for an arbitrary configuration of the random vector po- 
tential parameterized as aj(r) = Vj0(r). 

In a finite system, a degree of the wave functions' lo- 
calization (or a lack thereof) can be inferred from the 
inverse participation ratios 



V n =< 



J \%P(r)\ 2n dr 
L 2 (/|V(r')|W) 



L 2 



(23) 



where the Gaussian averaging over the disorder field 
</>(r) was performed with the weight P[0(r)] cx 
cxp(~Jdr(V 2 ^ 2 /2T v ). 
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The expression (23) is in a stark contrast with that in 
the short-range case where the zero-energy wave func- 
tions demonstrate a "pre- localized" behavior and the 
participation ratios exhibit a multifractal spectrum of 
anomalous dimensions V n oc L~ ari+hn — . 

Eq.(23) is suggestive of a strong localization of the 
zero-energy states, the apparent localization length being 
of order ~ r„ ' . Moreover, it is conceivable that all 

1 /2 

the states up to the energy of order ~ vT v remain local- 
ized. It should be noted, however, that by modeling RMF 
as a random variable with the variance (2) one misses 
out on the possibility of percolating "snake" trajectories, 
akin to those found in the case of non-relativistic fermons 
in a smoothly varying RMF with the correlation length 
greater than w/eia. 

As far as two-particle processes are concerned, a diffu- 
sive behavior of the corresponding amplitudes for |ex ± 
£2! rfa) can be studied in a more familiar framework of 
the non- linear tr-model (NLoM). In the metallic regime 
(cf 3> la), a derivation of the corresponding (supersym- 
metric) NLoM would closely follow the solution of the 
RMF problem for non-relativistic spinful fermions with 
a gyromagnetic ratio equal twr>!£. This model features 
the same unitary symmetry as that with only an orbital 
coupling to the long-range RMF—. 

In the opposite ("undoped") limit, e_p < j a r , a sys- 
tematic derivation of the NLoM is hindered by the lack 
of a suitable expansion parameter, for the bare conduc- 
tivity takes values of order unity (in units of ~ e 2 /h). 
For methodological purposes, however, one can still in- 
troduce a large number of valleys N, similar to the analy- 
sis of dirty <i-wave superconductors with tp = 0, thereby 
allowing for a formal 1/iV-expansion. 

This caveat notwithstanding, it was recently argued 
that in the case of short-range disorder the NLoM ob- 
tained by non-abelian bosonization could hold all the 



way down to the lowest energies^. It would be inter- 
esting to see if this conjecture can also be extended to a 
long-range RMF, in which case the corresponding local- 
ization scenario would fall into the so-called "CI" uni- 
versality class^. Such a behavior develops in the pres- 
ence of a non-vanishing inter-valley scattering (due to, 
e.g., a screened RP) in addition to the small-angle (intra- 
valley) one described by Eq.(2). It is, however, less likely 
that even in the absence of inter-valley scattering the 
system would show any scale-invariant behavior (univer- 
sality class AM), as it does in the short-range case^. 

In light of the above, the argument of Ref^ that 
RMF scattering suppresses any quantum coherence be- 
tween pairs of time-reversed trajectories (which can re- 
sult in either localizing or antilocalizing behavior, de- 
pending on the relative strength of the intra- vs inter- 
valley scattering^) would only apply in the case of a 
smoothly varying RMF. As the aforementioned NLoM 
analyses suggest, in the case of RMF described by Eq.(2) 
the onset of localization gets merely postponed until 
greater length scales {e Ka v/rfj instead of e Ka v/"f t J , pro- 
vided that the bare conductivity a 1). 

In summary, we study the behavior of two-dimensional 
Dirac fermions subject to long-range-correlated random 
scalar and vector potentials. In the ballistic regime of 
large quasiparticle energies, we obtain estimates of the 
scattering rates manifested by the Drude conductivity, 
dHvA/SdH and Friedel oscillations. The distinct energy 
(density) dependencies of such rates provide a means 
of ascertaining the dominant mechanism of elastic scat- 
tering in graphene. In the complementary low-energy 
regime, we find a signature of strong localization in a 
long-range RMF and discuss pertinent localization sce- 
narios. 

This research was supported by NSF under Grant 
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